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■ Abstract 

We compute the light-cone momentum distributions of pions in the nucleon in chiral effective 
theory using both pseudovector and pseudoscalar pion-nucleon couplings. For the pseudovector 
coupling we identify 5-function contributions associated with end-point singularities arising from 
the pion-nucleon rainbow diagrams, as well as from pion tadpole diagrams which are not present in 
the pseudoscalar model. Gauge invariance is demonstrated, to all orders in the pion mass, with the 
inclusion of Weinberg- Tomozawa couplings involving operator insertions at the ttNN vertex. The 
results pave the way for phenomenological applications of pion cloud models that are manifestly 
consistent with the chiral symmetry properties of QCD. 
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I. INTRODUCTION 



Contributions from the meson cloud of the nucleon to deep-inelastic scattering (DIS) were 

h ri 

first discussed in the early 1970s by Drell, Levy and Yan [1[ and Sullivan 2]. It was of little 
interest, however, until in the early 1980s a possible enhancement of the pion cloud of a 
bound nucleon was suggested as an explanation of the nuclear EMC effect {3, 4]. While that 
issue is still of interest, the modern importance of the pion cloud in high-energy scattering 
originates in the work of Thomas [5] , who showed (within a particular chiral quark model) 
that the pion cloud contribution provided a natural explanation for the observed excess of 
non-strange over strange sea quarks in the nucleon. This calculation also predicted an excess 
of d over u quarks in the nucleon sea, which was confirmed by the violation of the Gottfried 



sum rule observed in the NMC experiment at CERN 



6], and the d/u ratio measured in the 



Drell- Yan reaction in pd and pp scattering by the E866/NuSea collaboration at Fermilab 

Although this excess of d over u, as well as the s — s asymmetry [gj], was motivated by the 
physics of spontaneous chiral symmetry breaking in QCD, the early calculations were based 
on models whose connections to QCD were not manifest. It was later realized {9), however, 
that because of their origin in chiral loops, these contributions had a nonanalytic dependence 
on quark mass that could only be generated by a Goldstone mechanism, which placed these 
effects on a more rigorous theoretical footing. In fact, it is a model- independent consequence 
of spontaneous chiral symmetry breaking in QCD that d — u and s — s are nonzero. The 
only open question is how large these asymmetries actually are. 

This issue has become even more important in the last few years because of the widespread 
interest in the so-called five-quark components of baryon wave functions [lo| . In particular, 
there is some evidence that the Roper resonance may not be a simple excitation of a single 



valence quark in the nucleon but rather involves a large five-quark contribution [ll|-|l5| . This 
naturally leads to questions about such components in the wave function of the nucleon, 
and there have been suggestions that these may be sizable. Even though the pentaquark 
now appears to be defunct, legitimate questions about five-quark components in baryon 
spectroscopy remain, especially within models such as the chiral quark soliton model 16 1. 

In the context of low energy tests of the Standard Model this issue is also very topical, with 
a potential asymmetry between the s and s parton distributions potentially yielding a large 
correction to the value of sin 8\y derived from the NuTeV measurement [17J, Q. Model- 
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dependent estimates of the five-quark component of the nucleon wave function invo 
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vmg 



three light quarks plus ss pairs were recently used to investigate this correction 

In all of these examples, the contribution from the meson cloud of the nucleon has been 
a crucial factor in reconciling the physics results. Indeed, because of the model-independent 
constraints imposed through chiral symmetry on the nonanalytic behavior 9], the presence 
of meson cloud contributions to the various observables is firmly established, justifying 
the extensive theoretical attention that has been paid to this issue over the past decade. 
A further motivation relates to the fact that lattice QCD is now providing considerable 
information on the moments of parton distribution functions 23J , but at larger quark masses 
than occur in nature. The extrapolation of those moments as a function quark mass to the 



24, and 



physical point is critically dependent on knowing the correct nonanalytic behavior 
that in turn is uniquely determined by the pion cloud of the nucleon. 

Given the phenomenological importance of the meson cloud, it is unfortunate that the 
literature contains a number of sometimes contradictory results for meson cloud contribu- 
tions to DIS. In particular, the majority of calculations, from the original work of Drell et 
al. [l| and Sullivan |2], to the early studies of the nonanalytic behavior (9], and essentially 
all model analyses of pion cloud effects in DIS in between [251 ]. have used pseudoscalar (PS) 
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2811 . The restoration 



coupling, which is, by itself, inconsistent with chiral symmetry 
of chiral symmetry can be achieved through the addition of a scalar "cr" field; however, 
determining its practical consequences for DIS is problematic because of uncertainties in 
identifying the nature of low-mass scalars in meson spectroscopy [291 ]. 

More recently, the matrix elements of nonsinglet twist-2 operators were computed by 



Chen and Ji 



30 



3jj and Arndt and Savage [32J using lowest order, heavy baryon chiral 



perturbation theory. By construction this theory uses pseudovector (PV) tiN couplings, 



which are manifestly invariant under chiral transformations [28|, |33|, [34| . The twist-2 matrix 
elements are related through the operator product expansion to moments of parton distri- 
butions measured in inclusive DIS. At lowest order in the pion field 7r, the effective chiral 
Lagrangian for the interaction of pions and nucleons, consistent with chiral symmetry, can 



be written as 



28 



33 



34] 



C 7 



N 



1 



1p N ^ T ■ (7T X 7T) ijl. 



V: 



(1) 



where ipN is the nucleon field, /„- = 93 MeV is the pion decay constant, and qa = 1.267 is the 
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nucleoli axial vector charge. The first term in the Lagrangian ([[]) gives rise to the well-known 
"rainbow" diagram in which a pion is emitted and reabsorbed by the nucleon at different 
space-time points. The second term generates the so-called pion tadpole contribution, in 
which two pion fields couple to the nucleon at the same point. The tadpole term does not 
appear in the PS theory, which leads to some of the differences between the moments of 
twist- 2 parton distributions from pion loops computed in the PS [9] and PV 30-32] theories. 

In this paper we present a detailed analysis of the light-cone momentum distributions of 
pions, and the corresponding recoil baryons, inside a physical nucleon within the chirally 
symmetric effective field theory defined by the Lagrangian (CQ). We further contrast the 
results with those obtained in the PS theory, conventionally defined by the Lagrangian 

Cl% = -g^NN 4>n ij 5 r-7rijj N , (2) 

where the 7rNN coupling constant g n NN is related to qa and f n by the Goldberger-Treiman 
relation {35) ]. 

9-kNN 9a /„\ 
with M the nucleon mass. Unlike the earlier chiral effective theory calculations which only 

n 

computed the light-cone distributions of pions [3JJ or considered the nonanalytic behavior 



of their moments to lowest order in the pion mass [30|, 132J, we compute the complete set 
of diagrams relevant for DIS from nucleons dressed by pions resulting from the Lagrangians 
([1]) and without taking the heavy baryon limit. In particular, we demonstrate explicitly 
the consistency of the computed distribution functions with electromagnetic gauge invari- 
ance. This requires consideration of the Weinberg- Tomozawa terms 28|, |36], which although 
entering at higher orders in m n , are nonetheless essential for ensuring conservation of charge 
to all orders in m n . 

Formally, the light-cone distribution functions can be defined by considering the vertex 
renormalization constant Z-f 1 — 1 for the physical 7*iViV vertex. If the contributions from 
the pion cloud are not large, then Z\ ~ 1 and one has Z^ 1 — 1 fa 1 — Z\. The light-cone 
distributions fi(y) associated with a particular contribution % can then be defined as j^] 

{l-Z l ) i = I dyfi{y), (4) 







where y = k + /p + is the momentum fraction of the physical nucleon carried by the pion, with 
k and p the four-momenta of the pion and physical nucleon, respectively. Note that in this 



work we define the light-front "+" and "— " components of a four-vector as = v° ±v z . 
It will be most natural to perform the calculation in light-front coordinates; however, as 
demonstrated in Ref. 38[ for the nucleon self-energy, the results for the model- independent, 
long-distance physics associated with the pion cloud are reproduced in any formalism, be it 
instant form, covariant, or light-front. 

Because the tadpole diagrams in the PV theory involve pions emitted and absorbed at 
the same point, these can only contribute at zero light-cone momentum fraction, y = 0. 
A novel feature of the PV theory, however, is the appearance of ^-function contributions 
arising also from the rainbow diagram at y = 0. Inclusion of these singular terms is in fact 
vital for preserving gauge invariance in the PV theory. In the following we will elucidate the 
origin of the 5-function terms and discuss their possible physical implications. 

In the rest of this paper, we derive in Sec. [Til the light-cone momentum distribution 
functions for pions in a physical nucleon, including both the pion rainbow and tadpole 
contributions. The analogous nucleon momentum distributions are discussed in Sec. IHH 
where the Weinberg- Tomozawa couplings required by gauge invariance are introduced. The 
complete set of pion and nucleon distributions allows us for the first time to explicitly verify 
the symmetry relations between them. The application of the results to inclusive DIS is 
described in Sec. [IV] within the convolution approach, where we discuss the effect of the 
pion loop corrections, including the 5-function contributions, on various sum rules. Finally, 
we summarize our findings and their implications in Sec. [V] 



II. PION LIGHT-CONE MOMENTUM DISTRIBUTIONS 

Assuming isospin symmetry, the contributions from pion loops to nucleon matrix elements 
cancel for isoscalar combinations. In this work we therefore consider only isovector light- 
cone momentum distributions, fi %v \y) = ff{y) — f^iu)- The distributions are evaluated by 
taking the "+" components of the operators. We begin by considering the contributions 
from diagrams involving a direct coupling to the pion, as illustrated in Fig. [TJ 

The isovector light-cone momentum distribution function for the pion rainbow diagram 
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FIG. 1: Contributions to the pion light-cone momentum distributions in the nucleon, from (a) the 
pion rainbow, and (b) the pion tadpole diagrams. 



in Fig. Ufa) is given by 
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(5) 



where the pion and nucleon propagators are given by 



D n = k 2 -ml + ie, D N = (p — k) 2 — M 2 + ie, 



(6) 



respectively. In Eq. (jSJ) and throughout this work the spinors u(p) are normalized such that 
u(p) u(p) = 1. (Note that for the pion coupling diagrams, the isovector combination p — n 

" [3l|.) Using the Dirac equation, the distribution 



is equivalent to the notation u n + 



7T 



function in Eq. (jHJ) can be decomposed into several terms, 



1/ 



9A 



2fJ 7(27r) 



d 4 k 



AM 2 p ■ k 2M 2 p ■ k 

+ h 

DID 



N 



2y5[y 



yd 



(7) 



To evaluate the expression in Eq. ([7j), we proceed to first integrate over the k~ component 
of the pion momentum. Without loss of generality, we work in a frame where p± = 0. For 
the first term in (JJJ) we use Cauchy's integral formula and take the single nucleon pole by 
closing the contour in the upper half-plane, 

k\ + M 2 + is 



D 



N 



{p + - k A 



P 



hr 



0. 



(8) 



p + — k + 

Note that the arc contribution at infinity vanishes for this term. A straightforward calcu- 
lation then yields the contribution which we associate with the on-shell part of the nucleon 
propagator, 

9\M 2 r 2 y(k 2 ± + y 2 M 2 ^ 



f ion \y) 



; dk]_ 



[kl + y 2 M 2 + (l-y)ml] 



2 ■ 



(9) 



where for convenience we have factored out all isospin factors. This term corresponds to the 
usual pion distribution function for the "Sullivan process" [2j, which is computed from a PS 
ttN coupling. Performing the k± integration with an ultraviolet cut-off A gives 



g\M 2 



f {on \y) = -iT^y 



77-4(1 -y) fml(l -y) + M 2 y 2 ^ 1 f 1 r2 



+ lo 



772,2 (1 -y) + M 2 y 2 V A2 



A 2 " 



(10) 



where only the leading term in A has been kept. 

In the second and third terms in Eq. (J7|), the cancellation of the nucleon propagators 
leaves tadpole-like contributions involving pion propagators only. This will have important 
consequences for the structure and interpretation of the distribution functions. Integrating 
the second term in (J7J) proportional to y/D 2 over k~ gives a result which is zero everywhere 
except at k + = yp + = 0, 
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After multiplication by y, this term therefore vanishes. 

Evaluation of the third term in Eq. (j7]) proportional to (yp ■ k)/D 2 requires particular 
care. Since p-k — (p + k~ +p~k + )/2, integration over k~ involves terms such as in (TTT1) which 
vanish, as well as terms of the type J dk~ (k~ /D 2 ). To evaluate the latter, one can make 
use of the identity 

k~ d 1 



Dl dk+D w 



(12) 



to write 



dk-^f = J dk~jj- = 2tT7 log ) S(k + ), (13) 



where fi is a mass parameter. This gives rise to a 5-function contribution (again defined 
without isospin factors) 

f {5) (y) = ttA^ I dk\ log ( k l±^\ 5(y), (14) 



with the total pion distribution function then given by 



fW{v) = ±f [on) (y) + 4/ w (y)- (is) 
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CI 

f^ S \y) = 7 — A f N 2 m ^l°g m ^ $(y) + terms involving A, (16) 



When integrated over k± from to A, the function f^(y) becomes 

A 

4(47T/ ff ) 2 ' 

where we have isolated the leading nonanalytic term in m 2 , whose coefficient is independent 
of the short-distance physics parametrized by A. 

Thus the pion rainbow diagram has contributions from the regular y > region as well 
as the ultrasoft y = point. In the PS formulation of pion-nucleon interactions, where the 



pion light-cone momentum distribution is given by just the first term in 
term in ( fi~5l) does not arise, and the total PS pion distribution is given by 



5q. (J7D, the /W 

3 



In addition to the pion rainbow diagram in Fig. [H(a), there is also a contribution from 
the term in the Lagrangian (CQ) involving two pion fields, which gives rise to the tadpole 
diagram in Fig. [TJ(b). A straightforward calculation of the resulting contribution to the pion 
light-cone momentum distribution gives 

/&(!/) = Jl J^^P)(-^)<P)^ 2k+S (k + -yP + ) (18) 
i f d^k p ■ k ( k+\ 

A similar calculation to that for the f^(y) term above gives the final pion tadpole contri- 
bution as 

/&W = -J^JJ, J dki log (^^) S(y) = 2/^)(y), (20) 

where for convenience we have defined the distribution f( tad \y) without the isospin factor 
2. Integrating over k± from to A then yields a result for f^ tad '(y) similar to that for the 
fW{y) term in Eq. rtl4| ). 



/ (tad) (y) = 2(4?!/ y m K logm * S ( y } + terms involving A (21) 
2 



-j/^G/)- (22) 



;(*), 

Finally, combining the results in this section, the total isovector pion distribution in the 
nucleon, including both the pion rainbow and tadpole diagrams in Fig. [TJ can be written as 



r(iv ) ( , _ 4^M 2 [ m*(l-y) 

(47T/,) 



ml(l -y) + M 2 y 2 



+ log (ml(l - y) + M 2 y 2 ) 



+ , v 2 m^logm^ <5(y) + terms involving A, (23) 



where we have explicitly separated the long-distance contributions from the short- distance 
effects involving the cut-off A. In particular, the total term proportional to the ^-function is 
found to vanish in the limit g& — > 1, as one would have in a PS pion-nucleon theory 
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Note that a similar behavior for the S(y) term is observed for the pion cloud contribution 



to the isovector form factor 



39, 4o]- 



III. NUCLEON LIGHT-CONE DISTRIBUTIONS 



For the nucleon rainbow diagram in Fig. [2]^a), the distribution function for the neutron 
is 2 x that of the proton, so that for the isovector distribution one has (y) = ~f P (y)- 
Specifically, the isovector nucleon distribution function is given by 



fH v \y) 



-M 
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d A k 



u(p) (#75) 



D N 
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N 



(24) 



which, using the Dirac equation, can be written as 

2 , Mi 
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AM 2 (k - 2yp ■ k) AM 2 y 



D W D 2 N 



1 

EL 



D n D N 

The first term (~ 1/D^D 2 ^) in Eq. (1251) corresponds to the on-shell nucleon contribution, 
and is proportional to the function f^ on \y) in Eq. (Q. The second term (~ 1/D V D^) arises 
from the off-shell components of the nucleon propagator, while the third term involving the 
single pion propagator (~ l/D^) contributes only at y — 0, and is proportional to f^(y). 
The total isovector nucleon light-cone distribution function arising from the nucleon rainbow 
diagram is then given by 



fP(y) 



/M(y)-/W(y)+/W(y) 



[•(off) 



(S) t 



(26) 



where 



f {oS) (y) 



g\M 2 

(47T/ )r ) S 



dk\ 



y 



k\ + y 2 M 2 



y)m\ 



9\M 2 ( ml(l - y) + M 2 y 2 \ f M 2 
WfJ* yl0g { A 5 ) + ° U 



(27) 
(28) 



with the k± integration taken up to the cut-off scale A. 
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FIG. 2: Contributions to the nucleon light-cone momentum distributions in the nucleon, from 
(a) the nucleon rainbow, and (b) the tadpole diagrams, as well as (c) the contributions from the 
Weinberg- Tomozawa coupling required by gauge invariance. 

For the pseudoscalar model the nucleon light-cone momentum distribution is again just 
given by the on-shell term in Eq. ( l26l) , 

= -f (on) (y)- (29) 

This result was also obtained in the infinite momentum frame calculation of Drell, Levy and 
Yan 

On the other hand, the PV theory contains, in addition to the off-shell and ^-function 
pieces, also the contribution from the operator insertion at the NNtttt vertex, Fig. [2(b). 
The distribution function associated with this diagram can be written 

ftU(y) = -| J^u{p) 1 +u{ P )^-5{k+-y P +) (30) 
i [ d*k 1 / k+\ 

= -e 7 (So 1 isT _ fJ- (31) 

Performing the k~ integration and again using the relation in Eq. (fl3"|) then gives 

/#id)(v) = -2/ (tad) (2/)- (32) 

Comparing the two tadpole contributions, Eqs. (|20|) and (|32l) . one finds that their sum 
vanishes, 

fSUto) + fmUv) = o- (33) 
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Tadpoles themselves therefore have no net effect on the sum of the pion and nucleon light- 
cone distributions, and hence on the total nucleon charge, as will be discussed in Sec. HVl 
below. 

Finally, the light-cone momentum distribution associated with the Weinberg- Tomozawa 
diagrams in Fig. Etc), which arises from the derivative coupling in the PV theory, is given 
by 

2 , ,4, 



/$(y) = 4M(^f 



9a 



2fJ J (2tt) 
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(35) 



After k integration, one then obtains 



4%) = 4/w(y) - %f w (y)- 



(■(off) 



(36) 



One observes then that the sum of the pion (Eq. ([IB]) ) and WT (Eq. fl3"6"|) ) distributions is 
proportional to the nucleon light-cone distribution (Eq. (|2"6"|) ). 



/f)(y) + /™(y) = 4/W(y) + 4/W(y) - 4/W(y) 



e(on) 



(off). 



-(S)l 



*fP(y), 



(37) 



which, as will be discussed in the next section, is a necessary condition for the gauge invari- 
ance of the theory. 



IV. SUM RULES 



The effect of pion loops on observables such as parton distribution functions can be 
computed by considering DIS from a nucleon viewed in the infinite momentum frame 
Here the total contribution to the nucleon isovector quark distribution qffl from the pionic 
corrections in Figs.[T]and[2]can be written in convolution form (for a review see e.g. Ref. j^]), 

qP(x) = Z 2 + (if ) ® i V) ) (*). (38) 

i 

where gjy is the isovector quark distribution in the bare nucleon state No, and the 
sum over i includes the nucleon and pion rainbow, tadpole and WT contributions. The 



11 



symbol ® denotes convolution; for the pion rainbow contribution, for instance, one has 
(f^ ®<}k ){x) = $l{dy/y) f^ V \y) Q^ V \^/y), where is the isovector quark distribution 
in the pion, and similar for the other terms. While the shapes of the bare quark distribution 
functions, and especially for the WT and nucleon tadpole terms, are not a priori known, 
charge conservation is nevertheless ensured as long as the distributions are normalized ac- 
cording to 

dxqf v \x) = 1, (39) 



o 



for all i = N, 7T, WT, nN. 

The wave function renormalization factor Z 2 in Eq. ( 138|) can be evaluated from the nucleon 



self-energy £ by 



4l|-|44| 



*r'-i = £ 



(40) 

p =M 



where the derivative is taken at the nucleon pole, po = M. For weak pion fields, explicit 
evaluation gives 

Z -i _ i „ i _ z 2 = 3 f X dy (y) + f (oS) {y) _ /W(y) j j (41) 

where the factor 3 arises from the sum of the proton and neutron intermediate state con- 
tributions, for either a proton or neutron initial state. In fact, Eq. (14 ip is required by the 
Ward-Takahashi identity which relates the wave function and vertex renormalization factors 
by Z 2 = 3(1 — Z±), where Z± is the lowest moment of the nucleon distribution function /at 
for a proton initial state. 

Integrating the convolution expressions in Eq. (I3"8j) over all x and using the relation in 
Eq. ( 1371) one obtains, for a nonsinglet quark distribution q^ v \x) normalized to unity, the 
quark number sum rule, 

\xqf\x) = i = z 2 + 3 / dy {f {on) (y) + f ioS \y)-f {5 \y)), (42) 



o 

with the contributions from the /^ad) anc ^ ^i(tad) terms canceling as in Eq. (I3"3"|) . From the 
expression for Z 2 in Eq. (14ip one can verify that the baryon number of the nucleon is not 
modified by pion loops. 

Inclusion of the y = contributions is vital for formal sum rules to be satisfied; if these 
were to be evaluated for y > only, the 6(y) terms would be missed. Indeed, as was 
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illustrated in Ref. 



45] in the case of the Bukhardt-Cottingham sum rule for the twist-3 



polarized parton distribution Q2{x), if such S(x) terms are present then sum rules based on 
the formal operator product expansion appear to be violated when evaluated for x > only. 
The same was observed for the corresponding sum rule for the chirally odd twist-3 parton 
distribution h^x). In the context of dispersion relations, 5(x) terms reflect the presence of 
subtractions. A deep-inelastic scatterin experiment to test such sum rules would thus find 
them to be violated, since x = is never measured. In the convolution calculation the S(y) 
terms would thus imply a 8(x) contribution to the quark distribution qf^\x), which would 
lead to the violation of both the Gottfried |46| and Adler sum rules j^J , if these are applied 
to nonzero x only. 

The appearance of the 5(x) terms in the isovector quark distributions is possibly an 
artifact of taking chiral perturbation theory too literally in this context. For example, 
although a tadpole diagram inevitably leads to an s-independent contribution to the forward 
Compton amplitude — and hence a subtraction in the corresponding dispersion relation - 
it may no longer be justified to assume a point-like ttN four-point coupling that gives rise to 
the tadpole when s 3> m n . Nevertheless, even if the Adler sum rule is in the end satisfied, 
the presence of 8{x) terms in chiral perturbation theory indicates that in order to confirm 
this sum rule it is essential to include the small- x regime. If chiral perturbation theory is 
taken literally, a 8(x) contribution will presumably soften into a peak near x = in QCD. 
Of course, none of these considerations apply to the case of the PS pion-nucleon coupling, 
which contains nonzero contributions only at y > 0. 



V. CONCLUSION 



In this paper we have presented a detailed derivation of the isovector light-cone momen- 
tum distributions of pions in the nucleon for both the pseudovector and pseudoscalar ttN 
interactions. In the PV theory the direct coupling to the pion gives rise to a pion rainbow 
distribution, fi tv \ and a distribution from the pion tadpole, /f^d)' wr th the latter involving 
emission and absorption of a pion at the same point. Consequently the tadpole contribu- 
tion to the momentum distribution exists only at zero light-cone momentum, y = 0. We 
found, however, that the pion rainbow diagram also has a 5- function component, on top of 
the regular terms at y > 0. The S(y) contributions arise because the PV ttNN vertex is 
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proportional to the momentum of the pion, which leads to a more singular integrand in the 
k~ — > oo limit than in the case of the momentum-independent PS coupling. 

In addition to the pion distribution, we have also computed the isovector light-cone 
momentum distribution of the corresponding recoil nucleon dressed by the pion, including 
contributions from the rainbow diagram with nucleon coupling, f^"\ and a tadpole coupling 
at the tttcNN vertex, f^n^, as well as the Weinberg- Tomozawa terms J^t which arise at 
the same order in the chiral expansion of the PV Lagrangian. The latter are in fact essential 
for preserving the gauge invariance of the PV theory. The nucleon rainbow distribution 
contains pieces that are associated with the off-shell components of the nucleon propagator, 
in addition to the on-shell components, and a ^-function term analogous to that in the pion 
light-cone distribution /i . In contrast, the nucleon distribution in the PS model only 
contains the on-shell contributions. The WT distribution fyfy has contributions from the 
off-shell components of the nucleon and from the singular <5-function terms. The singular 
nucleon tadpole distribution is related to the pion tadpole by /^ad) = — /jv(tad) ^ or an V- 

Combined, the sum of the pion rainbow and WT distributions in the PV theory is shown 
to be proportional to the nucleon rainbow distribution, 

fi iv) (y) + f&%) = -4/#°(v), (43) 

which is a necessary condition for ensuring gauge invariance. Within the convolution model, 
we have demonstrated that this in fact gives the correct number sum rule: the sum of all the 
pion cloud contributions does not alter the valence quark number of the nucleon, although 
it does of course affect the shape. This is also true in the PS formulation. 

Our results pave the way for future phenomenological applications of pion cloud models 
that are manifestly consistent with the chiral symmetry properties of QCD. Most calculations 
in the literature of the effects of the pion cloud in deep-inelastic scattering have been made 

|25[, for which chiral symmetry is not manifest without 
the addition of a scalar field. In fact, the absence of the 5-function contributions in the 
PS version of the rainbow diagrams leads to the coefficient of the leading nonanalytic term 
in the chiral expansion of the quark distribution function moments that differs by a factor 



using the PS tcN coupling 



4/3 from the PV result [37|]. A detailed study of the phenomenological consequences of the 
distributions derived here in the PV theory will be the subject of an upcoming study j^ . 



14 



Acknowledgments 



This work was supported by the DOE contract No. DE-AC05-06OR23177, under which 
Jefferson Science Associates, LLC operates Jefferson Lab, DOE contract No. DE-FG02- 
03ER41260, the DOE Science Undergraduate Laboratory Internship (SULI) Program, and 
the Australian Research Council through an Australian Laureate Fellowship and the ARC 
Centre of Excellence in Particle Physics at the Terascale. 



[1] S. D. Drell, D. J. Levy and T. M. Yan, Phys. Rev. D 1, 1035 (1970). 

[2] J. D. Sullivan, Phys. Rev. D 5, 1732 (1972). 

[3] M. Ericson and A. W. Thomas, Phys. Lett. B 128, 112 (1983). 

[4] C. H. Llewellyn Smith, Phys. Lett. B 128, 107 (1983). 

[5] A. W. Thomas, Phys. Lett. B 126, 97 (1983). 

[6] P. Amaudraz et al, Phys. Rev. Lett. 66, 2712 (1991); M. Arneodo et al, Phys. Rev. D 50, 

Rl (1994); M. Arneodo et al, Phys. Lett. B 364, 107 (1995). 

[7] E. A. Hawker et al, Phys. Rev. Lett. 80, 3715 (1998). 

[8] A. I. Signal and A. W. Thomas, Phys. Lett. B 191, 205 (1987). 

[9] A. W. Thomas, W. Melnitchouk and F. M. Steffens, Phys. Rev. Lett. 85, 2892 (2000). 

[10] R. Vogt, Prog. Part. Nucl. Phys. 45, S105 (2000). 

[11] O. Krehl, C. Hanhart, S. Krewald and J. Speth, Phys. Rev. C 62, 025207 (2000). 

[12] W. Melnitchouk et al, Phys. Rev. D 67, 114506 (2003). 

[13] D. Guadagnoli, M. Papinutto and S. Simula, Phys. Lett. B 604, 74 (2004). 

[14] N. Mathur et al, Phys. Lett. B 605, 137 (2005). 

[15] M. S. Mahbub et al, Phys. Lett. B 679, 418 (2009). 

[16] K. Goeke, M. V. Polyakov and M. Praszalowicz, Acta Phys. Polon. B 42, 61 (2011). 

[17] D. Mason et al, Phys. Rev. Lett. 99, 192001 (2007). 

[18] W. Bentz, I. C. Cloet, J. T. Londergan and A. W. Thomas, Phys. Lett. B 693, 462 (2010). 

[19] J. Alwall and G. Ingelman, Phys. Rev. D 70, 111505 (2004). 

[20] M. Wakamatsu, Phys. Rev. D 71, 057504 (2005). 

[21] Y. Ding, R.-G. Xu and B.-Q. Ma, Phys. Lett. B 607, 101 (2005). 



15 



[22] W.-C. Chang and J.-C. Peng, Phys. Rev. Lett. 106, 252002 (2011). 

[23] G. S. Bali et al, Phys. Rev. D 86, 054504 (2012). 

[24] W. Detmold, W. Melnitchouk, J. W. Negele, D. B. Renner and A. W. Thomas, Phys. Rev. 

Lett. 87, 172001 (2001). 

[25] J. Speth and A. W. Thomas, Adv. Nucl. Phys. 24, 83 (1997). 

[26] M. Gell-Mann and M. Levy, Nuovo Cim. 16, 705 (1960). 

[27] M. Gell-Mann, Phys. Rev. 125, 1067 (1962). 

[28] S. Weinberg, Phys. Rev. Lett. 18, 188 (1967). 

[29] M. R. Pennington, Mod. Phys. Lett. A 22, 1439 (2007). 

[30] J.-W. Chen and X. Ji, Phys. Lett. B 523, 107 (2001). 

[31] J.-W. Chen and X. Ji, Phys. Rev. Lett. 87, 152002 (2001); Erratum-ibid. 88, 249901 (2002). 

[32] D. Arndt and M. J. Savage, Nucl. Phys. A 697, 429 (2002). 

[33] E. E. Jenkins and A. V. Manohar, Phys. Lett. B 255, 558 (1991). 

[34] V. Bernard, N. Kaiser and U. G. Meissner, Int. J. Mod. Phys. E 4, 193 (1995). 

[35] M. L. Goldberger and S. B. Treiman, Phys. Rev. 110, 1178 (1958). 

[36] Y. Tomozawa, Nuovo Cim. A 46, 707 (1966). 

[37] C.-R. Ji et al, in preparation (2012). 

[38] C.-R. Ji, W. Melnitchouk and A. W. Thomas, Phys. Rev. D 80, 054018 (2009). 

[39] M. Strikman and C. Weiss, Phys. Rev. C 82, 042201 (2010). 

[40] N. Kaiser, Phys. Rev. C 68, 025202 (2003). 

[41] J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields, Mcgraw-Hill (1965). 

[42] M. B. Hecht, M. Oettel, C. D. Roberts, S. M. Schmidt, P. C. Tandy and A. W. Thomas, Phys. 

Rev. C 65, 055204 (2002). 

[43] M. Alberg and G. A. Miller, Phys. Rev. Lett. 108, 172001 (2012). 

[44] C.-R. Ji, W. Melnitchouk and A. W. Thomas. larXiv:1206.367H [nucl-th]. 

[45] M. Burkardt and Y. Koike, Nucl. Phys. B 632, 311 (2002). 

[46] K. Gottfried, Phys. Rev. Lett. 18, 1174 (1967). 

[47] S. L. Adler, Phys. Rev. 143, 1144 (1965). 

[48] K. Hendricks et al, in preparation (2012). 



16 



